Extinction in Discrete, Competitive, 
Multi-species Patch Models 

David Chan, Matt McCombs, Sarah Boegner, Hye Jin Ban, 
Suzanne L. Robertson 

Department of Mathematics and Applied Mathematics, 
Virginia Commonwealth University, Richmond, Virginia 

November 1, 2013 
Abstract 

In this paper we extend the results of Pranke and Yakubu in [5] for 
extinction in discrete competitive patch models. For a system of n species 
on m patches, we define conditions under which one species is a "superior 

competitor" to another and show that this is sufficient for one species to 
drive another to extinction. We also illustrate the result with an example 
for three species on three patches. 

1 Introduction 

Dispersal between patches has been shown to play a significant role in the dy- 
namics of populations that arc distributed over multiple habitats or patches, 
often stabilizing complex dynamics [1,7,9]. Hastings [6] showed that passive 
dispersal will never be selected for in a temporally constant environment, even 
if there is spatial variation among patches. However, dispersal may be advan- 
tageous in a temporally varying environment [8]. 

Pranke and Yakubu analyzed discrete-time competitive systems in [2,3], in- 
cluding dispersal of one species in [4]. In [5], Frankc and Yakubu explored a 
discrete competitive system of n species dispersing between 2 patches. They 
found conditions under which one species could drive another to extinction. 
Here wc extend their results to n species on m patches and provide an example 
showing the importance of dispersal to the dynamics of the model. 

2 Preliminaries 

Consider the following system, 

(Xi(t + l),X2(i + 1), . . . ,Xm(i + 1)) = F(xi{i),X2(i), • . • ,X,x>(i))> (1) 
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for n species on m patches which takes the form of n ■ m equations. For each 
l<i<m, l<j<n, and (xi, X2, . . . , Xm) G 3?+ x 3?!Ji x • • • x 3fi!Ji (m-times) 
where Xi = {xi^\,Xi^2, Xi^„) we define Xij{t+1) = Fj_j(xi, X2, . . . , Xm) where 
Fij : Ji^ X X • • • X Ji^ In particular Fij takes the form 



where Xij{t) is species j in patch i at generation t, gij : 3?" SR+ are monotone 
decreasing functions. Q j.j.fc is the competition coefficient in patch i for species 
k with respect to species j, and 6j € [0, 1] measures the diffusion of species 
j between the patches. Sj = 0 corresponds to no diflFusion and when Sj = 1, 
species j disperses so population sizes arc equal across all patches. We define 
i^* as the map F composed with itself t times. 

Consider for species j, 6j = 0 where there is no diffusion between patches. 
This would simplify _F!ij(xi,X2, . • • ,Xm) = ^i.j(xi) = .T,;j9ij(xi). Further if we 
assume that there were no other species on patch i then we have Fij{xij) = 
Xijgij{xij). This is a one-dimensional map. Let x* be the unique positive fixed 
point for this map which exists by assumption. If 0 < a; < x*, then Fij{x) > x. 
If a; > X* , then Fij{x) < x. So under iteration of Fij, lij = Fij{[0,x*]) is the 
largest attracting invariant interval in 3?+. The following shows that based on 
these competitive assumptions on the Qij the populations will not grow without 
bound. 

Lemma 1: In system (2), no point has an unbounded orbit. 

Proof: Let (xi,X2, . . . ,Xm) € 9?" x 3?!f: x • • • x 3?" . To show that orbits are 
bounded we need to show that the sequences {i^j*^(xi,X2, . . . jXm)}^^^ are 
bounded. To do this we first show that for each species j, 





2 



By definition 

-Pij(Xl,X2,...,Xm) 



/ n ^ 



~ ^^'i ( XI <^h,3,kXh,k {t) 



\k=l 



If 



Kk=l 



then 



-Fi,j(xi,X2,...,Xm) < Xij{t)gij ^Y^^'^ij^kXi,k{t)j 

= max < Xij{t)gij I ^ aij^kXi,k{t) 



If 



then 



-Fi,j(xi,X2,...,Xm) < Taaxxisj{t)g0j [^ai=ij^kXi3,k{t) 
= max I Xij{t)gij I ^ ai^j,kXi,k{t) 



\.k=l 



Thus niaxj{Fj,j(xi,X2, . . • ,Xm)} < maxj{a;ij(t)3ij (Efc=i 
Further because each gtj is monotone decreasing 

with ttij j = 1. Recall that if Xij < maxlij, then Fij{xij) < max/ij . Also 
if Xij > max/ij, then Fij{xij) < Xij. Thus for each l<i<m, l<j<n, 
and each point (xi,X2, . . . e 3?" x x • • • x we have max{Fjj} < 
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maxij {xij,maxlij}. Therefore the sequences {F*^(xi,X2, . . . ,Xm)}j>^ are 
bounded for 1 < i < m and 1 < j < n. This completes the proof. 



Additionally we have the following: 

Lemma 2: In system (1), the limit set of every point in 3?" x x • • • x 5i" 
is a nonempty subset of the compact invariant set B D W, where B = 
{(xi,X2, . . . ,Xi„) e 3?" X 5R!J: X • • • X 5f!J:| for each 1 < j < n, maxj{xij} < 
maxj{max Jjj}}, and W = {(xi, X2, . . . , Xm) G 3?" x 3?" x • • • x 5R!J:| for each 
1 < i,k < m and 1 < j < n, m-<m-i)S - ^hji^ + 1) < ^kjit + 1) < 



TO-(m-l)5, 

m— (m 



Proof: Let (xi, X2, • • • , Xm) e 3?" x 5R" x • • • x 3?" then from Lemma 1 we 
have that a;(xi,X2, . . . jXm) e 3?" € B. To show that ti;(xi,X2, . . . ,Xm) € 
W^, we observe the following where Sj € (0, 1). For convenience we will write 
9i,j {J2k=i 0ii,3,kXi^k{t)) as Qi^j (•). Without loss of generality consider species j 
on patch 1 we have 

Fj,j(xi,X2,...,Xm) = ^1 - Xij{t)gij{-) 

S ' 6 ' 

+ ^a;2,j(t)ff2j(-) H h ■^Xm,j{t)gm,j{-) 

Next multiplying by ^™ -'^^'^j yields 

m - (m - l)dj ^ \ - (m - (to - 1).^,)^ , , 

-fi,jiXi,X2, . . . ,Xn,j — Y Xij[t)gij[-) 

3 3 

m — (to — 1)5,- (5, 
+ ^^-^;^-2.(tK.(-)H--- 

m — {m— l)5j 6j , . 

+ ■■■H ~^rn,j\t)gm,j{') 

Since > 1 we have I^^ilpiliifi > ii and "^-^7-^^^^- (l - ^1^6,) = 

(m-(m-i)6j) ^ ^-j^ _ rn^^^y comparing term by term with F^j^ we have 



5jm 

that Ffej < as well as Fkj > ^_^^_^s, PiJ tlius 

F <p TO - (to - 

TO-(TO-l)^/''^- 5,- 
Finally since Xij{t + 1) = j(xi, X2, . . . , Xm) we obtain 

6j m — {m— l)Sj 

~Xi^j ^ Xk,j ^ 7 Xi^j- 



m — {m— l)5j 5 
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This completes the proof. 

3 Results 

For this model we assume that each species will diffuse between the m patches 
at a particular rate, given by the Sk's. For species k, if the rate Sk = 0 then 
there is no diffusion, where if ^fe = 1 then the populations will diffuse so that 
the populations in each patch will be equal. We show that in this model a 
species can drive another species to extinction if particular conditions are met. 
We define these conditions here. 

For system (1) and 6k,Sn G (0, 1], a species k is superior to species n if on 
some patch i and for some constant c > 0 for all points (xi, X2, . . . , Xm) € BdW, 
we have 



We note here that these conditions depend on the linear combination of the 
growth rates and well as the diffusion rates. Also note that for a diffusion rate 
of 1 these conditions simplify to having the average growth rate of the species 
k being larger that the average growth rate of species n. 

Theorem 1: For system (1), if species k is superior to species n, then species n 
will go extinct. 

Proof: Note that by Lemma 2 we have that the u;— limit set of every point is a 
nonempty subset of the compact invariant set B fl W. Let (xi,X2, . . . ,Xm) G 
B nW. We will show that (xi , X2 , . . . , Xm) 0 as t,; x for 1 < i < m. 

Suppose inequality (i) is satisfied on patch i so species k is superior to 
species n. Consider the function V : 5i" x 3?^ x • • • x 5f" 5f+ where 
y(xi,X2, . . . ,Xni) = We show that V decreases along the orbit so that 

xi,n = 0 for 1 < Z < m, thus V is a Lyapunov function. 




or 
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Let xi^k, xi^n > 0 for 1 < / < m, then we have 

y(f(xi,X2,...,X^)) ^ (1 - 5r,)Xi,n{t)gi,n (•) + ^ ELl Xh,n{i)9h,n (•) 

F(xi,X2, . . . ,x^) ((1 _ 5k)Xi,k{t)9i,k (•) + ^ Er=i Xh,k{t)9h,k (•))' 

Using the inequality rn-{m-i)5j ^hj(^ + 1) ^ + 1) < ""^T,"^^*^' + 1) 

we get 

y(-F(xi,X2,...,Xm)) 
F(xi,X2,...,Xni) 

(1 - 6n)Xi,nm,n (•) + ^ ( "'-^r'^'" ) ^r^l (■) 

((1 - ^fe)ara(i)5a (•) + ^ ( ^-(^'■-1),, ) Er=i XiMt)9h,k (•))' "^^'^ 

(1 - ^n)ffi.n (•) + ^ ( "-^r '^'" ) ^r=l 9H,n (Q 
((1 - (5fe)5»,fe (•) + '± { m-(^-l)S, ) ELl 9h,k (•)) 

So from inequaUty (i) we have that ^^^^'^^''^^^ ' ' ' ' •'^m)) ^ ]^ 

t/(xi,X2, . . . ,Xm) 

On the other hand, suppose inequality (ii) is satisfied on patch i so species 
k is superior to species n. Consider the function V : x 5R" x • • • x 5i!J: — ^ 5i+ 

where y(xi,X2, . . . ,Xm) = ''" ■ Let xi^k,xi,n > 0 for 1 < Z < m, we 

have 

V(f (xi,X2,...,Xm) ^ Ellll XLn{t)gi,n (') (E" 1 ■^^^fc)'' 

y(xi,x2,...,xm) (Ellia^^i.fcWs!,^ E]lia;j,n 

Again using the inequality ar/^jft+l) < a;;j(t+l) < Xij{t+ 

1) we get 

y(F(xi,X2,...,x^) Xi,n9iA-) + ( "'"^Z"^^'^" ) El^i,;^»a;i,„ft)g;,„ (■) 



nxi,X2,...,x^) + ( „,_(^-_,),J E"M^.^.feWg^fe (•))' 

+ 5n ) Ei = l,i^i (•) 



(9i,k{-) + (^rn-im-l)6k) ^1=1,1^^1 9l,k (O) 

m 

m— (m— l)(5n 
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From inequality (ii) we have that w^"'^"'"' ' -^t")) ^ ^ Thus from either 

l/(xi,X2, . . . ,Xm) 

inequahty we obtain that F is a decreasing along orbits for (xi, X2, . . . , Xm) € 
BnW wit'hxi^k,xi,n > Ofor 1 < / < m. Let (pi,P2, • • • ,Pm) e cj(xi,X2, . . • ,Xm), 
then by definition there exist a sequence <i — )■ 00 where F**(xi,X2, . . . ,x„i) — >■ 
((Pi) P2j • • • , Pm)- We want to show that F**^(xi, X2, . . . , x^) — *• 0 as — )• 00. 

If Pi k, Pin > 0 then we have that w^^"*"' " " ' ^"l^^ < 1. This cannot 

^(Pl,P2,---,Pm) 

hold for an w— limit point, which implies that Pi^n = 0 and Pi^k > 0. If Pi^n > 0 
and pi_k = 0 then the sequence {-F*' (pi, P2, ■ • ■ , Pm)} would be unbounded 
which is a contradiction. From Lemma 2 we conclude that if Pj^„ = 0 then 

Pl^n = 0 for 1 < Z < TO. Thus w(xi,X2, . . . ,Xm) C {(xi,X2,. . . , X^) € 

5R!f. X 3?'^ X • • • X ^l\xi,n = 0, 1 < ? < to}. 

Finally, consider the points that do not enter iJ n on which species k is 
superior to species n. Since each condition is a strict inequality we have that 
these hold on a neighborhood of BnW. From Lemma 1, we have that every point 
must enter this neighborhood that on every orbit, V is eventually decreasing 
and the result holds for every point with xi^k,xi^n > 0. This completes the 
proof. 



4 Example 

In this section we illustrate the importance of dispersal with a three-species, 
three-patch example. Wetake gi^- = exp ^r^j — X]fc=i "^jj-fc^i,*:) withru = 1.5, 
r22 = 1.4, — 1.3, and all other r^j = 1.1, i ^ j, and all aij — 0.1. If there is 
complete dispersal, so 5i = 62 = S3 = 1, then the average growth rate of species 
1 is the highest and it is a superior competitor to both species 2 and species 3, 
driving them to extinction by Theorem 1 (with c — 1). 

Whcm there is no dispersal, so 61 = S2 = 63 = 0, a different species persists 
in each patch (species 1 in patch 1, species 2 in patch 2, and species 3 in patch 
3) driving all other species to extinction. If we set 5i = 62 = S3 = S, then as 

5 increases from 0 to 1 we see coexistence for low S, sec Figure 1. Increasing 
further we observe species 3 becomes extinct and species 1 and 2 coexist on each 
patch. For S closer to 1 only species 1 persists on all three patches. It is clear 
through this example that dispersal has strong influence on the dynamics of the 
entire system. 
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Figure 1: Model (1) with gij = exp (^nj - J2l=i(^t-j,kXi.kj where rn = 1.5, 
= 1-4, = 1.3, — 1.1, i ^ j, and aU aij = 0.1 All species disperse at 
the same rate (Si = = 5;^ ^ S). Equilibrium values for each species plotted 
against 6. 
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